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We propose a novel representation of cosmic microwave anisotropy maps, where each multipole 
order £ is represented by £ unit vectors pointing in directions on the sky and an overall mag- 
nitude. These "multipole vectors and scalars" transform as vectors under rotations. Like the 
usual spherical harmonics, multipole vectors form an irreducible representation of the proper ro- 
tation group 5*0(3). However, they are related to the familiar spherical harmonic coefficients, 
aim, in a nonlinear way, and are therefore sensitive to different aspects of the CMB anisotropy. 
Nevertheless, it is straightforward to determine the multipole vectors for a given CMB map and 
we present an algorithm to compute them. A code implementing this algorithm is available at 
http://www.phys.cwru.edu/projects/mpvectors/ Using the WMAP full-sky maps, we perform 
several tests of the hypothesis that the CMB anisotropy is statistically isotropic and Gaussian ran- 
dom. We find that the result from comparing the oriented area of planes defined by these vectors 
between multipole pairs 2 < £i 7^ ^2 < 8 is inconsistent with the isotropic Gaussian hypothesis at 
the 99.4% level for the ILC map and at 98.9% level for the cleaned map of Tegmark et al. A partic- 
ular correlation is suggested between the ^ = 3 and £ = 8 multipoles, as well as several other pairs. 
This effect is entirely different from the now familiar planarity and alignment of the quadrupole 
and octupole: while the aforementioned is fairly unlikely, the multipole vectors indicate correlations 
not expected in Gaussian random skies that make them unusually likely. The result persists after 
accounting for pixel noise and after assuming a residual 10% dust contamination in the cleaned 
WMAP map. While the definitive analysis of these results will require more work, we hope that 
multipole vectors will become a valuable tool for various cosmological tests, in particular those of 
cosmic isotropy. 



I. INTRODUCTION 

A great deal of attention is currently being devoted to 
examining the power spectrum of the cosmic microwave 
background (CMB) temperature anisotropics AT{n)/T 
extracted from the Wilkinson Microwave Anisotropy 
Probe(WMAP) 0,13, 1^,4 and other CMB data HllE 
H, 13, 0, 0| • Decomposing the temperature in spherical 
harmonics 

^^^Y^ae^YU^) (1) 

£rrt 

and deducing the angular power spectrum 
1 ^ 

m——£ 

as a function of £ allows cosmologists to fit the pa- 
rameters of cosmological models to unprecedented ac- 
curacy, possibly even probing the physics of the in- 
flationary epoch. Similarly, the power spectrum of 



the temperature-polarization cross-correlation function 
is teaching us about the physics of the reionization of 
the universe, presumably by the first generation of stars; 
and the power spectrum of the temperature-galaxy cross- 
correlation function is teaching us about the statistical 
distribution of matter in the universe. But does the only 
physics lie in the angular power spectra? Is the sky sta- 
tistically isotropic, so that any variation in the values 
of the individual multipole moments aem with fixed £ 
represents only statistical fluctuations, or could there be 
subtle correlations between the a^m? If the sky is statis- 
tically isotropic, is it Gaussian — are the aim of each fixed 
£ drawn from a Gaussian distribution of variance that is 
only ^-dependent? Are there other interesting deviations 
from the simplest picture? 

In the standard inflationary cosmology the answer to 
the question just posed is that in the linear regime, i.e. 
at low £, the aim are realizations of Gaussian random 
variables of zero mean, with variances that depend only 
on £ (statistical isotropy). This paradigm is so strongly 
believed, both because of considerable observational ev- 
idence and considerable theoretical prejudice, that rela- 
tively little (though some, e.g. attention has been 
paid to searches for deviations from statistical isotropy. 



*More information, color pictures, and code are available at 



|http : //www ■ phy s . cwru.edu/projects/mpvectors/ 



In this paper we set out to search for one particular 
deviation — special directions on the sky. We do this by 
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first constructing from eacli multipole moment 

mm ^ 

m— — ^ 

of tlie CMB sky a set of ^ unit vectors {u^^'*) | i = 1,...,^} 
and a scalar ^'^^ tliat completely characterize that mul- 
tipole. We then examine the correlations between pairs 
of such sets of vectors, and {{)(^2,i2)|^ compar- 

ing them with Monte Carlo simulations of CMB skies 
with statistically isotropic Gaussian random aim ■ If the 
sky is statistically isotropic with Gaussian random aim, 
{aimai'm') = CeSu'Smm'- Since the {{)(^'''')} depend only 
on the ai-rn, {w'-^^'*^-'}, and {{;(^2,*2)| should be uncorre- 
cted for £i^e2- 

We constructed these vectors for 2 < £ < 8 for a set 
of full-sky maps including the WMAP Internal Linear 
Combination (ILC) JlJ, and the WMAP map as cleaned 
by Tegmark et al. [13 • We applied four statistical tests 
to the set of vectors from these full-sky maps. We find 
that one of the tests is inconsistent with the hypothesis 
of statistical isotropy and Gaussianity at the 99% con- 
fidence level. This work complements recent work by 
Eriksen et al. [T^ looking at north-south asymmetries in 
N-point functions, work by Park 15] looking at genus 
curves, work by Hajian et al. 16] using the so-called 
Ki test [T^ which finds violations of statistical isotropy 
for 15 < £ < 45, and work by Vielva et al. 18] using 
the spherical mexican hat wavelet technique where they 
found a strong signal for non-Gaussianity. 

II. MOTIVATION FOR A NEW TEST OF 
STATISTICAL ISOTROPY AND GAUSSIANITY 

Tests of non-Gaussianity, as opposed to statistical 
isotropy, have a long and rich history. Motivation for 
those tests originally came from the realization that 
non-Gaussianity is a signature of structure formation by 
topological defects ; while inflation predicts Gaussian 
CMB anisotropies. Subsequently it has been realized 
that, even if inflation seeded the structure in the universe, 
CMB non-Gaussianity may be present as a signature of 
features in the inflationary model [20ll2ll2a I2I |2^. 
Finally, late-time processes in the universe will induce 
non-Gaussianity on small scales [25II2II27II2II29II30I] . 
The tests of non-Gaussianity include studies of the bis- 
pectrum and skewness 112, IIE yj. Ill, |3a, 113 , trispec- 
trum [3M Isp. Minkowski functionals and the genus 
statistic |l5n40l M, M, ii] j_|pherical wavelets [H 
m Ei], a combination of these pl E l lillsollsH. and 
many other methods |U l54L l55l]. Of particular in- 
terest was the claim for non-Gaussianity in the COBE 
4-year data [56|, but this was shown to be an artifact of 
a particular known systematic [H^]- Nevertheless, efforts 
to test the Gaussianity of the CMB continue, and most, 
though not all (e.g. 0; Elj): have so far given results 
entirely in agreement with the Gaussian hypothesis. 



As these previous studies have shown, it is a chal- 
lenge to test such fundamental assumptions as statistical 
isotropy and Gaussianity without theoretical direction 
on what deviations to expect — they can be violated in 
a very large number of ways each of which could easily 
be hidden from the test that was actually performed. An 
instructive example is the effects on the CMB of any non- 
trivial topology of the universe, for example a universe 
which is a 3-torus. If the length scale of cosmic topol- 
ogy (for example the length of the smallest non-trivial 
closed curve) is sufficiently short, this will manifest it- 
self in various two-point temperature-temperature cor- 
relations, such as the so-called "circles-in-the-sky" sig- 
nature [53: if there are closed paths shorter than the 
diameter of the last scattering surface, then the last scat- 
tering surface will self-intersect along circles. These cir- 
cles can be viewed by an observer from both sides — from 
one side in one direction on the sky, and from the other 
side in some other direction. The temperature as a func- 
tion of location around the circle as seen from the two 
sides will be very strongly correlated. One can there- 
fore search for such pairs of circles. A definitive direct 
search is currently being conducted by the original pro- 
ponents of the signature. However, for us this example 
serves to show precisely why it is so difficult to perform a 
comprehensive test of statistical isotropy and Gaussian- 
ity. The topology-induced temperature correlations are 
strong only on or very near the matched circle pairs, thus 
the multipole coefficients aim, and other statistics that 
are weighted averages over the entire sky, would be poor 
tools for searching for such circles. Thus testing these 
phenomena is in part a matter of continually searching 
for (preferably) physically motivated ways in which they 
manifest. 

The CMB data itself may provide motivation for 
searching for deviations from the standard infiationary 
predictions (of statistical isotropy and Gaussianity), es- 
pecially at large angular scales. An absence of large an- 
gular scale correlations in the CMB sky relative to the 
inflationary prediction was first noted by COBE DMR 
in their first year data |6Cll ]. which showed what was re- 
ported as an anomalously low quadrupole, C2. Because 
the cosmic variance in the quadrupole is quite large, it 
was widely dismissed as a statistical anomaly. The re- 
sult persisted and was strengthened by the COBE DMR 
four year data [gJ. The recent WMAP analysis shows 
a marked absence of power on scales extending from 60° 
to 180° to an extent that cannot be explained solely by a 
low quadrupole Note that the estimator used by the 
WMAP team has been shown to be non-optimal when 
applied to incomplete maps of the sky. When an alter- 
native estimator is applied [Tsl l6^ or a full-sky map is 
analyzed 0] the discrepancy becomes less significant. 

An absence of power on large scales is expected in some 
topologically non-trivial universes. In a compact uni- 
verse there is a spectral cutoff of long wavelength modes 
leading to a suppression of power near this cutoff. One 
method of looking for such a cutoff is the "circles-in-the- 
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sky" signature as noted above. In general, such compact 
topologies would lead to "special directions" in the uni- 
verse. To search for special directions we need a method 
of defining our directions. Our definition, as discussed in 
the next section, is to decompose the multipole into 
i unit vectors, these vectors are then studied to search 
for peculiar alignments. These vectors contain the full 
information of the ae,m but encode it in a different way 
that allows for one to more easily look for special direc- 
tions. In particular, the components of these vectors are 
non-linear combinations of the aim (for 9- fixed t). 



III. DEFINING "SPECIAL DIRECTIONS" 

One attempt to look at the statistical isotropy of the 
CMB on large scales was the analysis of the quadrupole 
and octupole moments of the WMAP sky by de Oliveira- 
Costa et al. [fi^ . They found that the quadrupole was un- 
usually small, and that the octupole was unusually planar 
and unusually aligned with the quadrupole. They identi- 
fied an axis with each multipole by finding, for each the 
axis ill around which the angular momentum dispersion 



\ — 

(4) 

is maximized. (Here aimiw) are the spherical harmonic 
coefficients of the CMB map in a coordinate system with 
its z-axis in the hi direction.) They found that the n2 
and the 113 directions: 



na (-0.1145,-0.5265,0.8424) 
hs ^ (-0.2578,-0.4207,0.8698) 



and 



(5) 



are unusually aligned — their dot product is 0.9838. This 
has only a 1 in 62 chance of happening if n2 and 113 are 
uncorrelated and the dot product is uniformly distributed 
on the sky. De Oliveira-Costa et al. |^ point out that 
these values of n2 and 113 could be explained by a universe 
which has a compact direction parallel to 112 and n3 and 
of length approximately equal to the horizon radius; but 
this is ruled out by other tests, including the absence of 
matched circles in these directions. 

It has been noticed that the quadrupole and the oc- 
tupole in the cleaned WMAP skies remain dominated by 
a hot and a cold spot in the Galactic plane — one in the 
general direction of the Galactic center, and the other in 
the general direction of the molecular cloud in Taurus. 
This raises the possibility that the observed correlation 
is dominated by foreground contamination. One would 
like therefore to examine in more detail the correlations 
between the aim corresponding to possible preferred di- 
rections, or correlations in directions between the various 
multipoles. 

The question then is how best to associate directions 
with the CMB multipoles. De Oliveira-Costa et al. |6^ 



associated only one direction with each multipole, corre- 
sponding to two real degrees of freedom, whereas the aim 
of a given £ have 2£ -I- 1 real degrees of freedom. The 
multipole, fe{^), in the multipole expansion of a function 
/(ri) on a sphere 

f{n) = J2fii^) = J2 aerr^Yimin) (6) 

e im 

can be fully represented by a symmetric, traceless rank 
£ tensor, Fi^,,,i^ [ik — 1,2,3). Such a tensor can readily 
be constructed from the outer product of £ unit vectors, 
and a single scalar, A'^^y (Strictly speaking these 
are headless vectors, i.e. points on the projective two- 
sphere. The sign of each vector can always be absorbed 
by the scalar. The sign of the scalar takes on physical 
signficance when we define a convention for the multi- 
pole vectors, such as, all of them point in the northern 
hemisphere.) 



A. Vector Decomposition 

The correspondence between these vectors and the 
usual multipole coefhcients can readily be seen for a 
dipole. A dipole defines a direction in space — the line 
along which the dipole lies. The standard correspondence 
is 

Yifl^z, Yi^±i^T^{x±iy). (7) 

Thus 



^ ai 

,ra 

m— — 1 

= A(i)(w(1'1\{)(1'1),i)(1'1)) • (sin 6* cos 0, sin 6* sin ^, cos 6*) 

= -e, (8) 

where e is the radial unit vector in spherical coordinates. 
For a real valued function, the vector's components are 
found to be 

^^■i) = -N/2a'i';i, 4i'i)=V2a™, = ai,o, (9) 

and A'^^^ = (which can then be used to construct 

the unit vector, v'^'^'^'^). 

To extend this to the £^^ multipole, we want to write 
heuristically that 

^ a,mYimm « . e) . . . (^(«) • e) , (10) 

m=-l 

for each of the £ directions given by e. This cannot be 
quite right since the product of £ vectors would contain 
not only components of angular momentum £, but also 
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of angular momenta (. — 2, £ — 4, etc. However, a sim- 
ple power-counting shows that, once the reality condi- 
tions have been imposed on 0£m, l unit vectors and a 
scalar contain the same number of degrees of freedom as 
does aim, namely, {21 + 1) real degrees of freedom. We 
therefore expect that the components of lower angular 
momentum found in the right hand side of equation Hl()|l 
are not independent. We shall see this explicitly in equa- 
tions (|13|l and H15() below, and more elegantly in section 
IIII Bl For now, let us treat equation Hl()|l as motivation 
and proceed. 

Instead of solving H1U|) directly for all -O^^'*) we peel off 
one vector at a time, finding first a vector (with 
components v[['^''), and a rank l— \ symmetric, traceless 
tensor, a^^'^\ We can think of ii as running over x, ?/, 
and z, or more conveniently over —1, 0, and 1. Similarly, 
we can write a^^'^-* as a 3 x 3 x • • • x 3 ((f — 1) terms) ma- 
trix, alj'^.^if' however, this hides its traceless, symmetric 
nature and makes it appear that a'^'^^ has far more in- 
dependent degrees of freedom than it actually does. It is 
therefore more instructive to write the 21— \ independent 
components as a^^^^^, with m — —{£ — 1), {£ — 1). 

We repeat this procedure recursively on the remaining 
symmetric, traceless tensor from the previous step. Thus 
we next peel from a*^^'^-' a vector, v^^''^\ and a rank i — 2 
symmetric, traceless tensor a^^''^\ and repeat until we 
have found the full set of £ vectors, {u^^'*) | i = 1, ...,£}. 
The scalar A^^\ is found in the last step when the sec- 
ond to last vector is peeled off. In this case the 
remaining symmetric, traceless tensor, a*-^'^"^^ is rank 1, 
and is just the product of the final unit vector and 
the scalar A'-^^ 

To apply the recursive procedure outlined above we 
use the following recursion relation to peel off one vector 



and that 



for j 



1,0,1 where 



(11) 



^±1 



{£ - m){£ + m) 
47rY {21- l)(2/-t- 1) 



{e,m) 



D 



±1 




3 {£±m-l){£±m) 
8^\l {2l-l){2l + l) 

3 l{£-m-l){£ + m-l) 
i^V (2;-3)(2;-l) 

3 l{£^m-l){£Tm) 
8^' 



(12) 



(2;-3)(2Z-l) 

For a given £ we peel off a vector using 
e (f-i) 1 



e-2 



(13) 



(14) 



The second term on the right-hand side of equation H13|l . 
is required to guarantee that this rank £ — 1 tensor is 
traceless. In other words, it subtracts off the trace. We 
can see this is necessary from the y^_2,m term in the 
recursion relation Hll() . The presence of this £—2 term is 
as anticipated in the discussion following equation H1U|) 

Plugging in the recursion relation Ullfl yields 4£ — 1 
coupled quadratic equations that must be solved for the 



1 unknowns: v 



{(,1) 



a} 



we find 



1 



1 



i=-i 



1 



and b^^'P . From them 



{21 + 1 equations), 

{21-3), (15) 
(1 equation). 



These equations are easily solved numerically, though see 
appendix ^ 

Notice that, again as anticipated below equation H1U|I . 
the components, b^^, are not relevant for further calcu- 
lations; they are functions of the v^^'^^ and a^l^^\„_jj^ and 
are not independent. This means that the vectors v^'^ we 
calculate are indeed unique as claimed. 

The general procedure now follows: from the given 

i-im- continue using 
t-2 m- This is repeated until 



aim we construct and a^^'^^ 



atZn to find ^(^.2) and -^'^^'^ 

we find a'^^'^^ which gives the final two vectors ^(^'^"i) 

and The result of this -procedure for the CMB is 

shown in figure ^ 



B. Alternative derivation of vectors 

In practice we have implemented the procedure de- 
scribed above (appendix ^ and solved the set of equa- 
tions (|15|l for the analysis we have performed. A math- 
ematically more sophisticated decomposition procedure 
that leads to the same set of vectors without the need to 
calculate the bm' begins by recognizing that 

Mn) = F^%^^ {e^-&- ■ ■ ■ e^] ^ fI^..^P^---^ (16) 

Here e is again a radial unit vector, there is an im- 
plicit sum over repeated indices (ii, •••,?£), and the 
square brackets represent the symmetric trace free part 
of the outer product. For a general symmetric tensor. 
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FIG. 1: An image of the sky as decomposed into the £ — 2 through 8 muhipole moments based on the first year WMAP 
results Pi as cleaned by Tegmark et al. 13]. Shown are the 'Y^^aimXim and the vectors calculated for these multi- 
poles. The vectors are drawn as "sticks" since they only defined up to a sign (thus they are "headless" vectors). See 
|http: //www. phys . cwru. edu/projects/mpvectors/ for a full sized color picture. 



n2 L'?/2J 



(2^)! 



{2£-2ky. 



k=0 



k\{l - k)\{£ - 2k)\ 



X g{ili-2 , , , gi2k-li2k gi2k+l---it)pi---P2k (17) 



X 5pip2 ■ • • ^p2k-lP2k ' 

where there is an impUcit sum over the repeated indices 
Pi J ■■■,P2k, and denotes symmetrization of the enclosed 
indices : 



(18) 



(Ui is the group of permutations of the numbers 1, £.) 
The particular combination 

O'' = [e'^e'^ ■■ -e''] (19) 
simplifies considerably because e°-e''Sab = 1 



,N2 L^/2J 

E 

fe=0 



{2l-2k)\ 



k\{£ - k)\{£ ^ 2k)\ 



X ^(*l*2 . . . J«2):-l«2fcgj2fe+l . . . gij;) 



(20) 



More importantly, it is easily calculable because of the 
recursion relation: 



(21) 



Similarly, the individual -O^^'*), can be peeled off one by 
one by recursion: 



(^) 



.(1,1) (^,1) 

'^11 '"l2---Jf 



-(^,2) {£,2) 



(22) 



Jl,l) 



assuming that the Fll!.,^ can be calculated. But, these 
are easily represented as integrals over the sky (just as 
are the aim)'- 



(,) ^ (2^+l)(2^)! [ ATjn) , 



(23) 



Alternatively, there is an explicit (analytic) relation be- 

{£) 

tween the Yim and the O^^""^' , so that F-^',,^^ can be 

expressed in terms of aim- Thus, given F(^) calculated 
from the sky, equation H23|l is a sequence of 2^+1 coupled 
quadratic equations for the v^^'^^ (and A^^^) analogous to 
equation in which the have been eliminated. 



IV. SOURCES OF ERROR AND ACCURACY IN 
DETERMINING THE MULTIPOLE VECTORS 

With the actual full-sky CMB maps, such as those 
that we use, the main source of error is pixel noise which 
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accounts for the imperfections in measured temperature 
on the sky. Pixel noise depends on a variety of factors, 
one of which is the number of times a given patch has 
been observed. Pixel noise for WMAP is reported as 
Cpix = fo/ V A^obs , where ctq is noise per observation and 
iVobs is the number of observations per given pixel |64j . 
For the WMAP V-band map the reported noise per pixel 
is (To = 3.11 mK and the variation in the number of ob- 
servations of each pixel is moderately small. Neverthe- 
less, it is important to account for the inhomogeneous 
distribution of pixel noise, as we do in our full analysis 
in Sees . IVl and I VII f where we find that the inhomogeneity 
of the pixel noise does not significantly change the main 
results) . For the purposes of illustrating the accuracy in 
determining the multipole vectors, however, we assume 
a homogeneous noise with the measured mean value of 
-^obs — 490 observations per pixel. 

The calculation of the pixel noise is straightforward. 
For equal-area pixels (as in Healpix 65]) subtending a 
solid angle ripix and assuming a full-sky map we have 



(24) 



so that, using a^-^ = /^AT(17)^ \ « constant, we have 



{aima*i,m>) 



f^pix^fipixr,:„(i7)y,:„(J7) /(AT(i7))'\ 



f^pixCpix Sw 5mm' ■ 



(25) 



Using the V-band map parameters, the pixel noise for 

NSIDE=512 Healpix resolution is y^A~(Tpix = 2.7 x 10"'' 
mK. We adopt this quantity as an estimate of pixel noise 
for all maps we use. 



A. Accuracy in Determining the Multipole Vectors 

We would like to find how accurately the multipole 
vectors are determined. To do this, we add a Gaussian- 
distributed noise with standard deviation Uaf,^ to each 



(26) 



where Af{^, <j^) denotes a Gaussian random variate with 
mean /i and variance a^. We repeat this many times in 
order to determine the distribution in the noise-added 
multipole vectors and their spherical coordinates (0,0). 
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FIG. 2; The accuracy in 6 (filled circles) and (p (empty 
squares) of a chosen multipole vector as a function of noise 
added to the aem- Both the mean value and scatter in the 
shift of the angles are shown. The i-axis value of unity corre- 
sponds to the scatter in the airn due to WMAP V-band pixel 



The mean and standard deviations of 9 and as a 
function of (Jae^ (which is in units of WMAP V-band 



pixel noise y fipix fpix) are shown in Fig. |21 The two 
panels show the effects on the £ — 2 and £ = 8 vec- 
tors. Note that both bias and scatter in 9 and 4> can 
be read off this figure. It is clear that the vectors are 
not extremely sensitive to the accuracy in the aim, and 
are determined to about ±1 degree for noise which is of 
order pixel noise. If the noise is much larger, however 

(faf„ ^ lOy^Apix cTpix), the accuracy in multipole vec- 
tors deteriorates to the point that they are probably not 
useful as a representation of the CMB anisotropy. Ac- 
curate determination of 9 and 4> is expected to be espe- 
cially important for higher multipoles, where the number 
of vectors is large. 
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V. TESTS OF NON-GAUSSIANITY WITH 
MULTIPOLE VECTORS 

Now that we have developed a formahsm to compute 
the multipole vectors, we would like to test the WMAP 
map for any unusual features. Clearly, there doesn't exist 
a test that simply checks for the "weirdness" of any par- 
ticular representation of the map, in this case the multi- 
pole vectors. We can only test for features that we specify 
in advance. Here our general goal is to test the statis- 
tical isotropy of the map and search for any preferred 
directions. 

Our motivation to devise and apply tests of statistical 
isotropy comes from findings that the quadrupole and oc- 
tupole moments of WMAP maps lie in the plane of the 
Galaxy as discussed previously. One could in principle 
extend this and ask whether any higher multipoles lie 
preferentially in this (or any other) plane, and devise a 
statistic to test for this alignment. Clearly, the number 
of tests one can devise is very large. Here we would like 
to be as general as possible, and choose tests that suit 
our multipole vector representation. We choose to con- 
sider the dot product of multipole vectors between two 
different multipoles as described below. 

A. Vector Product Statistics 

A dot product of two unit vectors is a natural measure 
of their closeness. One test we consider is dot products 
between all vectors from the multipole £i with those from 
the multipole £2- Since our vectors are really "sticks" (i.e. 
each vector is determined only up to a flip v —v), we 
always use the absolute value of a dot product. Further- 
more, motivated by the fact that the WMAP quadrupole 
and octupole are located in the same plane and that their 
axes of symmetry are only about 10° of each other [6^ . 
we also consider using the cross-products of multipole 
vectors of any given multipole: if the vectors of £1 and 
£2 lie in a preferred plane, their respective cross products 
are oriented near a common axis perpendicular to the 
plane. Dot products of these two cross products would 
then be near unity. 

We generalize this argument and make the following 
four choices for our statistic, which we shall call S. For 
any two multipoles £1 and ^2, we consider the following. 

1. Dot products of multipole vectors • t;(^2j)|^ 
where £1 ^ £2, u'^^^'*-' is the i^^ vector from the 
£1 multipole, and the j^^ vector from the 
£2 multipole. We call this statistic "vector- vector" . 
For a given £1 and £2, there are clearly M = £1^2 
distinct products. This statistic tests the orienta- 
tion of vectors. 

2. Dot products \v^'^''''> ■ (wf^^j) ^ {) (^2,fc))|/|^ («2 j) ^ 
V -ftrhere where v ^^i'*-' comes from the £1 mul- 
tipole and V and v (^^.fe) from the £2 multipole. 
We call this statistic "vector-cross" . For a given £1 



and £2 and j ^ k, there are M = £1^2(^2- l)/2 dis- 
tinct products. This statistic tests the orientation 
of a vector with a plane. 

3. Dot products |(-0(^i>*) x • (vC^^.k) ^ 

{)(^2,ni))|/(-|(^(^l,i) X 1(^(^2, fe) X {;(^2,m))|-) 

where and -O^^i'^) come from the £1 multi- 

pole and g^j^^ ^(f-2,m.) QQj^g from the £2 mul- 

tipole. We call this statistic "cross-cross". For a 
given £1 and £2 and i ^ j and and k ^ there are 
M = £i{£i - 1)4(4 - l)/4 distinct products. This 
statistic tests the orientation of planes. 

4. Dot products |(i)(^i>*) x u J)) • (w (^^''^^ xv^^^'™))] 
where v '^^^'^^ and v '^^^'i^ come from the £1 multipole 
and i} (^^.fe) ^nd v (^2,™) come from the £2 multipole. 
We call this statistic "oriented area". For a given 
£1 and £2 and i ^ j and and k m, there are 
M = £i{£i - 1)4(4 - l)/4 distinct products. This 
statistic tests the orientation of areas. Notice it is 
similar to the previous test but the cross products 
are unnormalized. 

B. Rank-Order Statistic 

Having computed the statistic in question, we would 
like to know the likelihood of this statistic given the hy- 
pothesis that the aim are statistically isotropic. The most 
straightforward way, and possibly the only reliable way, 
to do this is by comparing to Monte Carlo (MC) realiza- 
tions of the statistically isotropic and Gaussian random 
agm- To be explicit, we provide the algorithm for com- 
puting the rank-ordered likelihood, and also illustrate it 
in the flowchart Fig. |3| 

(a) For £1 and £2 fixed, a statistic S will produce M 

numbers (dot products). We will use S to test the 
hypothesis that the multipole vectors come from 
a map that exhibits statistical isotropy and Gaus- 
sianity. Calculate the M numbers for this statistic 
for the WMAP map. 

(b) To determine the expected distributions for these 

M numbers we begin by generating 100,000 Monte 
Carlo Gaussian, isotropic maps; in other words, we 
draw the coefficients aim by assuming {airnO*fjm') — 
C'lSu'Smm' ■ We are assuming statistical isotropy 
and Gaussianity, and this is the hypothesis that 
we are testing. We add a realization of inhomoge- 
neous pixel noise, consistent with WMAP's V-band 
noise, to each MC map. For each MC realization, 
we then compute the multipole vectors for multi- 
poles £1 and 4, and M dot products of the statistic 
S. Because the vectors of any particular realization 
do not have an identity (e.g. we do not know which 
one is the "fifth vector of ^ = 8 multipole") and 
neither do the dot products, we rank-order the M 
dot products. At the end we have M histograms of 
the products, each having 100,000 elements. 
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Pick statistic S 
(a) 



[2 



Pick £i ^ £2 
(a) 



Calculate the M values for the data 
(a) 



Rank order >Cwmap in 
these MC to get Ri^ 
(e) 



Test set of ranks 



MC to find the 
distribution for £ 
(d) 



Q Statistic 
(V C) 



Use MC to find 
the M distributions 
(b) 



Define L 
Calculate £wmap 
(c) 



MC Q distribution for 
final probability 
(VI) 



FIG. 3: A flowchart of the algorithm wo apply in order to extract the likelihood of the statistic, S. The lower case bold faced 
letters, such as (a), refer to the itemized points in section V B where more information can be found about the step. Also, the 
(V C) and (VI) refer to the sections in the paper where the details of these boxes can be found. 



(c) Wc would like to know the likelihood of the M prod- 

ucts computed from a WMAP map. To compute it 
we use a likelihood ratio test, which in the case of 
a single dot product of vectors would be the height 
of the histogram for the value of the statistic 5* 
corresponding to WMAP relative to the maximum 
height. Since we have M histograms, the likelihood 
trivially generalizes to 

f, _ TT iVj.yyMAP 

i-WMAP — — (.^'j 

where Nj^y/MAP is the ordinate of the j'*^ histogram 
corresponding to WMAP's j"' rank-ordered prod- 
uct, iVj^niax is the maximum value of j^^ histogram, 
and the product runs over all M histograms. 

(d) Now that we have the WMAP likelihood, we need 

to compare it to "typical" likelihoods produced by 
Monte Carlo realizations of the map. (An alter- 
native approach, computing the expected distribu- 
tion of the likelihood from first principles, would be 
much more difficult since one would have to explic- 
itly take into account the correlations between M 
products.) To do this, we generate another 50,000 
Gaussian random realizations of coefficients a^^m 



and m, aud compute the multipole vectors and 
the statistic S for each realization. 

(e) Wc rank-order the likelihood £wmap among the 

50,000 likelihoods from MC maps to obtain its rank 

(f ) Finally wc go to step (a) and repeat the whole proce- 

dure for all pairs of multipok;s (£1. £2) that we wish 
to test. Only when we have the complete set do we 
assign a probability. 

The rank Re^^e^ gives the probability that the statis- 
tic SwMAP is consistent with the test hypothesis. For 
example, if the likelihood of 5wmap is rank 900 out of 
1000, then there is a 10% probability that a Monte Carlo 
Gaussian random realization of the CMB sky will give 
higher likelihood, and 90% probability for a lower MC 
likelihood. We say that the rank order of this particu- 
lar likelihood is 0.9. If our CMB sky is indeed random 
Gaussian, we would expect the rank-orders of our statis- 
tics to be distributed between and 1, being neither too 
small nor too large. Conversely, if we computed the rank- 
orderings for three different multipoles for a particular 
test and obtained 0.98, 0.99 and 0.95, (or 0.01, 0.05 and 
0.02, say), we would suspect that this particular test is 
not consistent with the Gaussian random hypothesis. 
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C. How unusual are the ranks? 

We now consider how to quantify the confidence level 
for rejecting the hypothesis of statistical isotropy of the 
aim- Let us assume that we have computed the rank 
orders for N different pairs of multipoles, and obtained 
the rank-orderings of Xi for the i*'' one, where < Xi < 1. 
We consider the following parametric test. 

The test anticipates that the ranks will be unusually 
high. Let us first order the ranks Xi in descending or- 
der, so that Xi is the largest and x^ the smallest. We 
calculate the following statistic 

/■I ryi rVN-i 

Q{xi,...,xn)^ Nl dyi / dy2 . . ■ dyN- 

(28) 

If the ranks Re-^,12 were expected to be uniformly dis- 
tributed in [0,1], Q would be the probability that the 
highest rank is greater than xi, and the second-biggest 
rank greater than X2, ■ ■ ■ , and the smallest rank greater 
than xn ■ However we do not expect that the ranks from 
Gaussian random maps are uniformly distributed, and 
we treat Q merely as a statistic. We then ask that given 
the (possibly very small) value of Qwmap , what fraction 
of Gaussian random maps would give an even smaller Q? 
That number is our probability, and is computed in the 
next section. Note that, although apparently difficult or 
impossible to evaluate analytically, the right-hand side 
of Eq. (|28|l can trivially be computed using a recursion 
relation as shown in appendix IbI 



VI. RESULTS 

Table I shows the final ranks for the vector-vector, 
vector-cross, cross-cross, and oriented area tests for pairs 
(£1,^2)- These results correspond to the full-sky cleaned 
WMAP map from Tegmark et al. (results for their 
Wiener-filtered map from the same reference are essen- 
tially identical). The ranks from the WMAP ILC full-sky 
map are similar, and to make the presentation concise 
we show the ILC map oriented area ranks in Table I, but 
otherwise only quote final probabilities for the ILC ranks. 
As discussed earlier, a detailed analysis of cut-sky maps 
will be presented in an upcoming publication. 

For computational convenience we considered only the 
multipoles 2 < ^1,^2 < 8; we discuss the upper multipole 
limit in section IVI CI Note that, for the vector-cross test, 
(£1,^2) and (^2,-^i) products are distinct and both need 
to be considered. 



A. Vector-vector, vector-cross, and cross-cross 
ranks 

Table I shows that the vector-vector ranks are dis- 
tributed roughly as expected, nearly uniformly in the 



interval [0, 1]. The vector-cross ranks, however, are start- 
ing to shows hints of an interesting feature that will be 
more pronounced later in the cross-cross and oriented 
area tests: ranks that are unusually high: seven ranks 
out of 42 are greater than 0.9. The probability of this 
happening, however, is not statistically significant and 
may be purely accidental. 

The first big surprise comes from the cross-cross ranks: 
the (£1 = 3,^2 = 8) rank is 0.99988. This means that 
only 6 out of 50,000 MC generated maps had a higher 
likelihood than the WMAP map! In other words, the 
84 cross-cross dot products computed for this pair from 
WMAP lie unusually near the peaks of their respective 
histograms, which, recall, are built out of 100,000 prod- 
ucts from MC map realizations. The violation of statis- 
tical isotropy and/or Gaussianity therefore manifests it- 
self by a particular correlation between the vectors which 
makes the statistic S'wmap "unusually usual". Since we 
checked that the distribution of Monte-Carlo generated 
{£1 = 3, £2 = 8) cross-cross ranks is uniform in the in- 
terval [0, 1], it is easy to see that the probability of this 
rank being this high (or higher) is 6/50,000, or 0.012%. 
Admittedly, we checked 21 such cross-cross ranks, which 
raises the probability of find such a result to 0.26%. If 
we include all vector- vector, vector-cross, cross-cross and 
oriented area ranks, the probability rises to a still rather 
low 1.3%. 

Note that this effect is very different from the now- 
familiar orientation of the quadrupole and octupole axes; 
the quadrupole-octupole alignment is quite unlikely and 
results in multipole vector products which preferentially 
fall on the tails of their respective histograms. This is 
confirmed by the actual ranks for (2, 3) and (3, 2) multi- 
pole vector-cross products, and also the (2, 3) cross-cross 
and oriented area products, all of which are fairly low 
(< 0.05). This is easy to understand: in the cross-cross 
product case, for example, the fact that the multipole 
vectors lie mostly in the Galaxy plane implies that their 
cross products are roughly perpendicular to this plane. 
The dot products of those are then, by absolute value, 
very large, and hence unusual. What we are seeing here 
is that the {£1 = 3,^2 = 8) cross-cross products from the 
WMAP full-sky map are unusually usual. 



B. Oriented area ranks 

As mentioned above, one of the cross-cross ranks was 
extremely high. A much bigger surprise is found when we 
examine the oriented area ranks (see also Fig. 01): 2 (out 
of 21) are greater than 0.99, a total of 5 are greater than 
0.9, and a total of 8 are greater than 0.8! These ranks 
are clearly not distributed in the same way as those from 
a typical MC map. 

To examine the probability of the ranks being this high 
we perform the parametric test described in section IV Bl 
and in appendix^ we compute the statistic Q which, for 
a distribution of ranks expected to be uniform in [0, 1], 
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Ranks of Products of Multipole Vectors 




Vector- Vector 
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0.34475 


1 R 


1 9'^f^9 


(3, 


5) 


15 


0.21511 


30, 15 


0.36963 


0.78578 


30 


0.37187 


30 


0.67870 


ou 


n 7fi9j^/i 

U. ( 0Zo4: 


(3, 


6) 


18 


0.31507 


45, 18 


0.26683 


0.54146 


45 


0.75052 


45 


0.93546 




n f^9j^'^Q 

u.oziooy 


(3, 


7) 


21 


0.98772 


63, 21 


0.85874 


0.57072 


63 


0.55147 


63 


0.73650 


UO 


U.OOOZ4: 


(3, 


8) 


24 


0.76120 


84, 24 


0.98578 


0.60408 


84 


0.99988 


84 


0.99656 




n Q77f\f'\ 


(4, 


5) 


20 


0.41209 


40, 30 


0.28221 


0.84716 


60 


0.54035 


60 


0.52936 


DU 


u.ooyoo 


(4, 


6) 


24 


0.68840 


60, 36 


0.58372 


0.86140 


90 


0.74826 


90 


0.62266 


QD 


7'^7f\9 

U. 1 O 1 


(4, 


7) 


28 


0.85008 


84, 42 


0.51404 


0.95584 


126 


0.53715 


126 


0.88992 


1 9f^ 
±zu 








■^9 
oz 




119 48 






1 fiS 


n 84 "^74 


1 fi8 


n QQfinfi 

u.yyuuu 


iDo 




(5, 


6) 


30 


0.82148 


75, 60 
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0.66327 
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0.86242 


(5, 


7) 


35 


0.86884 


105, 70 


0.56542 


0.96116 


210 


0.59483 


210 


0.68920 


210 


0.98440 


(5, 


8) 


40 


0.83380 


140, 80 


0.96812 


0.34287 
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0.30403 


280 


0.24449 
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0.33959 


(6, 


7) 


42 


0.03742 


126, 105 


0.13831 


0.02203 


315 


0.20221 


315 


0.97286 


315 


0.78414 


(6, 


8) 


48 


0.92760 


168, 120 


0.91468 


0.78058 


420 


0.72850 


420 


0.62894 


420 


0.66873 


(7 


8) 


56 


0.03238 


196, 168 


0.04440 


0.02060 


588 


0.09552 


588 


0.25367 
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0.24871 



TABLE I; Ranks of the vector- vector, vector-cross, cross-cross, and oriented area statistics for multipoles 2 < ^1,^2 < 8 for the 
four tests we consider as appUed to the Tegmark et al. [l^ cleaned map. Also listed are the oriented area ranks for the ILC 
map. M is the number of products for a given statistic for each {£1,12) pair. In section fVlBI we perform a parametric test to 
compute the likelihood of oriented area ranks being this high. 



would be the probability of the largest one being at least 
as large as the largest actual WMAP rank and the sec- 
ond largest being at least as large as the second-largest 
actual rank, etc. This statistic, applied to ranks from the 
Tegmark et al. cleaned map, is 

Q = lO-(6.i2_i.i6) (29) 

where the error around the mean value is estimated by 
repeatedly adding the pixel noise of the map to the pure 
extracted a£m, as in Eq. (|26|l . and estimating the effect 
on the products of multipole vectors and their ranks. 

However, we have to be cautious not to over-interpret 
Q as the final probability — it has to be compared to 
Monte Carlo probabilities computed under the same con- 
ditions to determine its distribution. To this end, we 
generate 10,000 additional MC Gaussian random maps 
and compute Q for each. It turns out that only 107 of 
them produce Q lower than the WMAP value in (|29|1 . 
This is further illustrated in Fig. [S] which also shows 
the error bars on the WMAP Q. The probability of 
WMAP oriented area ranks being this high, according 
to the Q-test, is 1.07%, which corresponds to a 2.6-cr (or 
98.93%) evidence for the violation of statistical isotropy 



and/or Gaussianity. The connection of this (nearly) 3- 
a deviation to the nearly 3-ct deviation represented by 
the {£1 = 3, £2 = 8) cross-cross rank remains somewhat 
unclear. 



C. Further tests 

We have performed a few tests to explore the stabil- 
ity of the oriented area result. First, we have varied the 
multipole range from the fiducial 2 < £ < 8; the re- 
sults are shown in Table II. Increasing the lower limit, 
imin, leads to the final probability of 1% and 5.6% for 
^min — 3 and 4 respectively. Therefore, evidence for the 
violation of statistical isotropy weakens, but does so rel- 
atively slowly. This shows that our main result does not 
completely hinge on the quadrupolc and octupolc. Same 
can be said for the upper multipole limit, which gives 
strongest results for the violation of statistical isotropy 
with £aia.K = 8, but with decreasing ^max the result does 
not immediately go away. Finally, we have checked that 
increasing ^max to higher values, up to 12, does not pro- 
duce new ranks that are unusually high. The correla- 
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FIG. 4: Ranks of the oriented area statistics for the Tegmark 
fuU-sky map. The mean values correspond to the actual ex- 
tracted a^m, while the error bars were obtained by adding 
the pixel noise (the error bars are not necessarily symmetric 
around the no-error values). Note that an unusually large 
fraction of the ranks are high. Also note the low value of 
^1 = 2, £2 = 3 rank which is due to the alignment of the 
quadrupole and octupole which has been noted earlier [ij. 



O.lTdust- Although we are adding a significant contam- 
ination (the remaining dust is expected to contribute 
no more than a few percent to the rms CMB temper- 
ature jlj), the high oriented area ranks do not change 
much (they actually shghtly increase), and the Gaus- 
sian isotropic hypothesis is still ruled out at the 99.3% 
level. Clearly, dust contamination does not explain our 
results. Further, to mimic remaining foregrounds due to 
an imperfect cleaning of the map we tested adding a syn- 
thetic random Gaussian map which contributed 10% of 
the RMS temperature. We find that the oriented area 
statistic still disagrees with the Gaussian isotropic hy- 
pothesis at the 99.4% level. 

Finally, we have repeated the analysis with several 
other available full-sky CMB maps. As mentioned ear- 
lier, both maps analyzed by Tegmark et al. give the 
same probability for the oriented area statistic. For the 
WMAP ILC map we find similarly high ranks, giving an 
even smaller value for our statistic: Qilc = 2.44 x 10~^, 
and only 62 MC maps out of 10000 have a smaller value of 
Q; therefore, the high ranks in the ILC map are unlikely 
at the 99.38% level, corresponding to 2.7-a. 



VII. DISCUSSION 



Varying the Multipole Coverage 


^mi n 


QwMAP 


./(QmC < QwMAp) 


2 


7.61 X IQ-^ 


107/10000 


3 


3.13 X IQ-^ 


105/10000 


4 


3.12 X 10'* 


565/10000 




QwMAP 


/(QmC < QwMAp) 


8 


7.61 X 10"^ 


107/10000 


7 


3.72 X 10'^ 


394/10000 


6 


3.62 X IQ-^ 


2079/10000 



TABLE IL Final probabilities of the WMAP oriented area 
statistic as a function of multipole coverage /min < I < ^max. 
The fiducial case is 2 < £ < 8 and we have shown how the 
results change if the lower and upper bounds are changed. 
/{Qmc < QwMAp) is the fraction of MC random Gaussian 
maps that give a value of Q smaller than WMAP. 



tions are therefore most apparent in the multipole range 
2 < (^1,^2) < 8. 

We next check if the correlations can be explained by 
any remaining dust contamination. We use WMAP's V- 
band map of the identified thermal dust; this map was 
created by fitting to the template model from ReL [66| . 
We assume for a moment that 10% of the identified con- 
tamination by dust had not been accounted for, and we 
add it to the cleaned CMB map, i.e. Ttot = Icmb + 



Does the apparent violation of statistical isotropy or 
Gaussianity that we detected have a cosmological origin, 
or is it due to foregrounds or measurement error? The re- 
sults presented in this paper refer to full-sky maps, and 
it is known that there are two large cold and two hot 
spots in the Galaxy plane, and that any result that de- 
pends on structure in this plane is suspect. Nevertheless, 
it is far from obvious that the result is caused by the 
contamination in the map for the following reasons: 

• As Tegmark et al. argue, their cleaned map 
agrees very well with the ILC map on large scales, 
although the two were computed using different 
methods. The results we presented indicate a vio- 
lation of isotropy and/or Gaussianity at > 98.9% 
confidence using either map. Furthermore, instru- 
mental noise and beam uncertainties are completely 
subdominant on these scales. 

• The results come from an effect different than the 
quadrupole and octupole alignment: the latter is 
fairly unlikely, as discussed in Ref. jH3 1 while we see 
a particular correlation between the vectors that 
makes our statistic S unusually likely. 

• Perhaps most importantly, our results are mostly 
(though not completely) independent of the 
quadrupole and octupole, multipoles that might be 
suspect. For example, the second-highest oriented 
area rank is {£1 — 4,^2 = 8). Furthermore, Ta- 
ble IVI CI shows that if we only use the multipoles 
4 < ^1,^2 < 8, the oriented area statistic still rules 
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FIG. 5: The statistic Q for the oriented area statistic, computed from the Tegmark et al. cleaned WMAP map, is shown by 
vertical line. The shaded region around it corresponds to the uncertainty due to pixel noise, while the histogram shows the 
distribution of the statistic for MC-generated Gaussian maps. Only 1.07% of MC values of Q are smaller than the no-error 
value of the WMAP Q. The same fraction for the ILC map (not shown here) is 0.38%. 



out the Gaussian random hypothesis at the 94.4% 
level. 

At this time it is impossible to ascertain the origin of 
the additional correlation between the multipole vectors 
that we are seeing. One obvious way to find out more 
about their origin is to Monte Carlo generate maps that 
are non-Gaussian or violate the statistical isotropy, ac- 
cording to a chosen prescription, and see whether our 
statistics, S'wmap, agree with the S computed from MC 
maps. Of course, there are many different ways in which 
Gaussianity and/or isotropy can be broken, and there 
is no guarantee that we can find one that explains our 
results. 

Another possibility is to cut the galaxy (or other pos- 
sible contaminations) from the map prior to performing 
the vector decomposition. There are two approaches we 
can take: (1) use the cut-sky aim to compute the mul- 
tipole vectors and the statistics S and compare those to 
S computed from cut-sky Gaussian random maps, or (2) 
reconstruct the true full-sky a^m and compare with full- 
sky Gaussian random maps. The latter procedure is pre- 
ferred as one would like to work with the true multipole 
vectors of our universe, but reconstructing the full-sky 



map from the cut-sky information is a subtle problem 
that will introduce an additional source of error. Never- 
theless, the total error with a ~ 10° cut may still be small 
enough to allow using the vectors as a potent tool for 
finding any preferred directions in the universe. We are 
actively pursuing these approaches at the present time. 

Finally, Park pl recently tested the full-sky WMAP 
maps using the genus statistic and found evidence for the 
violation of Gaussianity at 2-3 a level, depending on the 
smoothing scale and the chosen aspect of the statistic. 
Furthermore Eriksen et al. HI find that the WMAP mul- 
tipoles with ^ < 35 have significantly less power in the 
northern hemisphere than in the southern hemisphere. 
It is possible that these results and the effects discussed 
in this paper have the same underlying cause, but this 
cannot be confirmed without further tests. 



VIII. CONCLUSIONS 

The traditional Y^ra expansion of the sky has many 
advantages. For one, each set of Yim of fixed £ form an 
irreducible representation of the rotation group in three 
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dimensions 0(3); for another, more than two centuries 
of effort have lead to a rich mathematical literature on 
the Yirn, their properties, and how to efficiently calcu- 
late them. The coefficients, aim of a Yim expansion of a 
function on the sphere are readily calculated as integrals 
over the sphere of the function times the y^*^. In this 
paper we have considered a different basis of an equiv- 
alent irreducible representations of the proper rotation 
group, O^^^ - for each £, the traceless symmetric prod- 
uct of £ copies of the unit vector of coordinates e(Jl). 
These are merely linear combinations of the Y^m, and 
so share many of the properties of them, albeit with a 
more sparse mathematical literature explicitly dedicated 
to their properties. In particular the coefficients Pj^^"^,,^^ 
of a (!?(£) expansion of a function on the sky are, like 
the spherical harmonic coefficients, calculable as integrals 
over the sphere of the function times C(£). 

We have expressed the F^J,,^^ as symmetric traceless 

products oil (headless) unit vectors {v^^'*)} and a scalar 
A^^\ The are highly non-linear functions of the 

aim- Thus, while in principle they encode the exact same 
information, they may make certain features of the data 
more self-evident. In particular we claim that these "mul- 
tipole vectors" are natural sets of directions to associate 
with each multipole of the sky. A code to calculate mul- 
tiple vectors from CMB skies is available on our website 
at http : //www .phys . cwru. edu/projects/mpvectors/ 
We have obtained the multipole vectors of the CMB 
sky as measured by WMAP, as well as the oriented areas 
defined by all pairs of such vectors (within a particular 
multipole). We have examined the hypothesis that the 
vectors of multipole £ are uncorrelated with the vectors 
of multipole £' for £ and £' up to 8. We have done this by 
comparing in turn the dot products of the vectors from £ 
with those from £' , the dots products of the vectors with 
the unit normals to the planes, the dot products of the 
unit normals to the planes with each other, and the dot 
products of the normals to the planes with each other. 
We have found that while there is nothing unusual about 
the distribution of dot products of the vectors with each 
other, the dot products of the normals to the planes with 
each other (and, to a lesser extent, the dot products of 
the unit normals to the planes with each other) are in- 
consistent with the standard assumptions of statistical 
isotropy and Gaussianity of the aim- To quantify this 
inconsistency we compared the distribution of these dot 
products with those from 50, 000 Monte Carlo simula- 
tions and found that they are inconsistent at the level of 
107 parts in 10000 for the Tegmark et al. cleaned full- 
sky map and 62 parts in 10000 for the ILC full-sky map. 
These results are robust to the inclusion of appropriate 
Poisson noise. The sensitivity to a Galactic cut will be 
explored in a future publication, but preliminary results 
suggest that the results persist within the error bars, but 



eventually decline in statistical significance as the uncer- 
tainties increase with increasing cuts. 
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APPENDIX A: VECTOR DECOMPOSITION 
EQUATIONS 

The vector decomposition equations we derived H15|l 
can be recast in a numerically more convenient form. The 
equations, as written, involve complex valued coefficients. 
Here we rewrite these equations in terms of their purely 
real components. To begin we note that the spherical 
harmonics satisfy -„(n) = (-l)"r/„(r2). Thus the 
decomposition coefficients of a real valued function, such 
as AT{fl)/T, satisfy a*_„ = (-l)'"af„. This shows 
that all the information about the function is encoded 
in the real part of aio (the imaginary part is identically 
zero) and the real and imaginary parts of aim for 1 < 
m < £- These are the 2£ +1 independent components we 
use in the vector decomposition. We thus only need to 
solve (221) ior < m < £. 

For notational convenience we drop the {£) superscript 
on ai-i^m-j, bm', and v. It should be understood that 
these quantities are associated with a particular multi- 
pole and step in the recursive decomposition procedure 
as outlined in section lTlI Al The correspondence between 
the dipole and Cartesian coordinate directions O allows 
us to identify v = {v^i,vq,vi) with standard coordinate 
axes via 

1 

V-l = — [V^+lVy) , 

vo = w^, (Al) 

Vi = --^(z)^- Wy). 

Finally the real and imaginary parts of the aim. are 

O'Tm = 2 ("^™ + '^im) and a"° = — {aim - Ofm) - 

(A2) 

Applying ljAl|l and (|A2|) to the multipole vector decom- 
position equations ((T^ gives the following equations 



14 



im im 



&0 = bl' 
\v\ 



1 



rc rc - I /rT/-»(^iO) / rc 



1 



.im _ r)('.m; im - 1 n^^'™J /„im f, 1 „ra - X 



1 



V2 



J_^(£,m) / rc - , im - 

^2 1 l"£-l,m-l"a: + "£-l,m-l": 



(A3) 



— n^'^'™) /'n'^c 7-, 4- n"n 7-, ^ 

J_r)(^, ™) /'„im - _ rc - \ 

l"<?-l,m-l''2; "^-l,m-l'^yj ' 



Note that the equations for the hm> are identical to those 
for aira with inserted in place of C'^^''"^ Here 

\ < m < I and 1 < m' < ^ — 2. These equations involve 
only real quantities and can thus be easily coded and 
solved. These are the equations we have implemented to 
find the multipole vectors. 



APPENDIX B: PROBABILITY OF RANK 
ORDERINGS 

Consider N numbers Xi, where < 2;^ < 1, and order 
them in descending order, so that x\ is the largest and 
xat the smallest. Let us then consider a set of variates yi 
uniformly distributed in the interval [0, 1], and also order 
them in descending order so that yi is the largest one and 
?/Ar the smallest. We ask: what is the probability that y\ 
is greater than x\ , and that 7/2 is greater than X2, ■ ■ ■ , 
and the yN is greater than xn- 

The probability that yi is in the interval [a;i, xi +dxi], 
Pi{xi)dxi, is 



Vi{xi)dxi = ( ^ Jxj dxi 



N 



(Bl) 



and the probability that yi is larger than xi, Pi{xi), is 
obviously 



Pi{xi) = / Vi{yi)dyi 



(B2) 



Given that yi is greater than xi, the probability that 
the 2/2 is in the interval [2:2, X2 + dx2] 

V2ix2\x,)dx2^{^'';') (B3) 

and the probability that the largest y is greater than xi 
and the second-largest greater than X2 is 

rl ryi 

P2ixi,X2) = / Vi{yi)dyi I 7^2 (2/2^1)^2/2 



N{N -I) 



f dyx r y^-^dy2.m 



We can continue this argument for all other yi and Xi , in 
descending order in Xi. The final probability, the joint 
probability of i"^ largest y being greater than Xi for all 
I, is given by 



Pn{xi,X2, ■ ■ ■ ,xn) = N\ dy 



dy2 



dy 



N- 



(B5) 

We would like to evaluate this integral. Even though 
the result will obviously be a polynomial in Xi, there is 
a total of 2^ terms and it is difficult to do the book- 
keeping. However, there is a simple recursion formula for 
this integral. Assume, more generally, that we want to 
compute 



'AT 



1 . 

dyi / 

1 Jx 



yi rVN 

dy2... / dyNy%. (B6) 

X2 Jxn-1 



One can then perform the innermost integral, and this 
leads to the recursion relation 



Im — 



1 



I y^N-i 



(B7) 



We are left with two (A'' — l)-tuple integrals. There- 
fore, starting from the A'^-dimensional integral, we can 
recursively bring it down all the way to A^ = 1, at which 
point it is an easy one-dimensional integral 



dyi 2/f 



1 



/3+ 1 



1-x^i 



(B8) 



for the required /3. Using the recursion relation (|B7|1 . 
together with ljB8|) . we numerically compute the proba- 
bihty in 
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